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ABSTRACT 

A balloon-borne spark chamber magnetic spectrometer has been used 

to measure separate spectra of positrons and negatrons in two flights 

during Summer, 1974. The flights reached atmospheric depths of 1.9 and 
-2 

1.5 g cm after slow ascents which enhance the statistical accuracy of 

the atmospheric secondary suotractlon. The total electron flux ir 
-2 -1 -1 -1 

about 0.3 m s sr MeV between 70 and 800 MeV, and increases toward 

lower energies. The positron spectrum decreases sharply toward lower 

-2 -1 -1 -1 

energies from a value of about 0.08 m s sr MeV at 650 MeV, and 
only upper limits are obtained fir positrons below 200 MeV. The impli- 
cations of these data are examined with regard to the problem of solar 
modulation. At energies above 180 MeV, the spherically symmetric Fokker- 
Planck equation with a diffusion coefficient proportional to particle 
rigidity provides reasonable fits to both che positron and total electron 
data. At energies below 180 MeV the data are consistent with a continu- 
ation of the same diffusion coefficient and local source of negatrons, 
or a change in the diffusion coefficient to a constant value. 
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I . INTRODUCTION 

Cosmic ray positrons and negatrons in the energy range from 15 
to 800 MeV are of astrophysical interest because of their ability to probe 
solar modulatior. This paper presents measurements of these particles 
made with a balloon-borne magnetic spectrometer. Precise measurements 
of these particles have, in the past, been hampered by several factors. 

The overwhelming background of cosmic-ray protons and the desirability 
of separation into e + and e has required complex instrumentation. 

The interpretation of fluxes measured with balloon-borne instrumentation 
has required an accurate knowledge of the contribution due to secondary 
electron production in the atmosphere. This knowledge can be obtained 
by combining a measurement of positron and negatron fluxes as a function 
of atmospheric depth with an appropriate theoretical model. 

The data presented in this paper represent a reasonable optimiza- 
tion of those parameters given the constraints of balloon experimenta- 
tion. Very large balloons (30 million and 50 million cubic feet) 
provided float measurements under a minimum of residual atmosphere. 

Slow ascents were used to provide accurate experimental growth curves. 
These data were then analyzed using a recent theoretical analysis of 
atmospheric secondary electron production and propagation to extract 
the contribution to the measured fluxes from sources external to the 
earth. 


II. INSTRUMENT AND BALLOON FLIGHTS 
The detector system used in these flights is very similar to 
that described by Daugherty (1974) ar.d Daugherty, Hartman and Schmidt 
(1975). Therefore, only a brief summary is presented here, with more 
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detailed description of only those features which differ from the ear- 
lier version of the Instrument. 

The apparatus, which Is shown In Figure 1, utilizes magnetic 
core digital spark chambers to determine the trajectories of charged 
particles through the Instrument, before and after they pass the gap 
of a permanent magnet. The direction and magnitude of the magnetic 
deflection provides the sign and rigidity of each particle traversing 
the detector. The spark chambers are triggered by a counter telescope 
consisting of a directional gas Cerenkov counter at.i two plastic scintil- 
lation counters. In addition, two anti-coincidence scintillation coun- 
ters are used to veto events In which a particle fails to pass through 
the magnet gap. 

The Cerenkov co-nter contains Freon-12 at 1 atm. pressure, which 

2 

provides an energy threshold of about 21 me for a particle of mass m. 

In principle, the gas Cerenkov counter Is insensitive to all protons 
and heavier nuclei with rigidities in the range of interest here (<1 GV) 
as well as all upward moving particles within the telescope aperture. 
Occasionally, however, a noise pulse from one of the Cerenkov photomul- 
tiplier tubes may occur in coincidence with the passage of a low-energy 
or upward moving particle through the telescope. As these random coin- 
cidences can yield particle trajectories which mimic good events, the 
reliability of electron identification has been enhanced by the use 
of two plastic scintillation counters D1 and D2 at depths of '''l and 3 
radiation lengths of lead. For rigidities below 450 MV they serve pri- 
marily as range counters for nuclei since, for example, a proton must 
have a rigidity of at least 350 MV to reach Dl, and at least 440 MV to 
reach D2. In addition, a proton wnich reaches Dl but not D2 must pro- 
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duce a signal in D1 of at least 5X minimum ionizing, and a proton which 
penetrates both D1 and D2 produces signals in D1 and D2 which are pre- 
dictable from the measured rigidity. The two coincidence counters SI 
and S2 are also pulse height analyzed, and are Important in eliminating 
nucleon background above about 200 MV. These data, along with the spark 
chamber data, are telemetered to the ground in a PCM format and recorded 
on magnetic tape. 

The experiment was flown on two balloons launched from Fort Churchill, 

Manitoba, on July 15 and August 3, 1974. The flights were planned to 

provide accurate positron and negatron flux measurements not only at the 

balloon float altitudes, but also for all atmospheric depths less than 
-2 

about 100 g cm . These ascent curves are essential for the estimation 
of primary and atmospheric secondary contributions to the observed flux- 
es. 

Unusually slow balloon ascents(l-2. 5m sec S were used in both 
flights to enhance the statistical accuracy of the ascent curves. Pri- 
mary electrons of low rigidity are observable from the Churchill area 
only between about 1800 CST and 0600 CST (Jokipil et al., 1968; Webber, 

1968; Israel and Vest, 1969; Daugherty et al., 1975). To avoid effects 
of the change from return albedo electrons (daytime) to primary electrons 
(night time) each flight was launched during the early evening. This 
allowed the slow ascent to bring the balloon to float altitude well 
before the increase in geomagnetic cutoff around 0600. The actual 
launch times were approximately 1830 CST July 15 and 2030 CST August 
3, with float altitudes being reached only after balloon sunrise, about 
0400 CST for both flights. The time at float altitude before the geo- 
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magnetic cutoff change was only about two hours, which was less than the 
planned time, but nevertheless satisfactory. 

-2 

The balloons reached float altitudes of 142,000 ft (1.9 g cm ) 

-2 

and 154,000 ft (1.5 g cm ). For rigidities greater than 200 MV (above 
the variable geomagnetic cutoff), 18 hours of float data were obtained 
for the two flights combined. 

III. DATA ANALYSIS 

Processing and analysis of the data was in many respects similar 
to that described by Daugherty (1974) and Daugherty et al. (1975); only 
those procedures which have changed substantially will be detailed here. 

The data for a single event consists of a variable-length list of 
set core (spark) addresses from the spark chamber, and the pulse-height 
data from the four scintillation counters. The spark chamber data were 
first processed by a pattern recognition computer program which classi- 
fies events as "good" (to be included in further processing), "bad" (dele- 
ted), or "undecided" (the program unable to make a dec^lon). The events 
classified as "good" are required to show a single t traversing 
the instrument within the nominal aperture with no obvious deflections 
except as expected in the magnet gap region. Spurious set cores are 
ignored if they do not appear to form a second track. "Bod" events are 
those which, for Instance, show either no complete tracks or more than 
one track through the instrument, or which appear to scatter from a pole 
face of the magnet. Approximately 95% of the 100,000 events recorded 
during the two flights were classified as either "good" or "bad" by the 
pattern recognition program. The remaining "undecided" events were 
examined using an interactive graphics unit where they could be deleted. 



- 6 - 

or edited and accepted, by the operator. Most of the events examined 
In this way were deleted because no complete track was discernible or 
because the track was obviously outside the Instrument aperture. The 
fraction of the graphically edited events which were accepted (^20%) 
showed evidence of short delta ray tracks or low spark efficiency In 
some region of the spark chamber. The events accepted graphically were 
tagged and merged with the original "good" events, which were then 
subjected to a least-squares fitting program to determine the particle 
trajectorl-s. The rigidity of each particle was calculated from Its 
magnetic deflection using a detailed map of t^e magnetic field, Including 
the fringe fields. Each event was examined and tagged for anomalies 
such as missing one cf the counters in the system or a deflection not 
fully compatible with that expected in the magnetic field. 

The angular deflection zero-reference and resolution were obtained 
from a histogram of deflection angles for those events which showed 
approximately minimum Ionization in all four analyzed counters. A 
correction of 2 mllliradlans was found to be necessary, end *».- -Wlec- 
tion angle resolution was found to be about 2 mllliradlans FWHM. The 
highest rigidity used in this analysis (800 MV) produces an angular 
deflection of about 9 mllliradlans. 

In order to eliminate possible background from protons or other 
cosmic-ray nuclei, restrictions were placed upon the pulse-heights in 
the four analyzed counters. Since the energy deposited by a nucleus in 
each scintillator is uniquely determined if the rigiditv is known, the 
restrictions were initially formulated so that any events were assumed 
to result from electrons which did not conform to the pulse-heights 
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expected for a nucleus with a rigidity determined from the spark chamber 
analysis, with suitable allowance for both rigidity and pulse-height 
resolution. With these rather loose restrictions It was found that, 
for positrons with energies above 120 MeV, a significant background 
existed on SI vs. S2 pulse-height plots, which was not present in the 
corresponding negatron plots. Additional restrictions have thereofre 
been included which require the pulse-heights In the coincidence scin- 
tillators to be less than some (energy dependent) maximum value. This 
reduces the electron recognition efficiency (determined from the nega- 
trons) to about 882 above 120 MeV, compared with 97-992 at lov»r energies. 

The event rates for electrons in various rigidity Intervals were 
examined to find the time of the change of geomagnetic cutoff mentioned 
above. The change was apparent in all rigidity intervals below 200 MeV. 
The time of the change appears to be rigidity dependent, but is the same 
for positrons and negatrons in a given rigidity range. Positron and nega- 
tron fluxes were therefore combined as shown In Figure 2. for each flight. 
Although in both flights the cutoff change appears to occur substantially 
later at the higher rigidities, only the first two float intervals of 
each flight were used for rigidities below 200 M e V. The entire float 
period was used for rigidities above 200 MeV. 

The observed positron and negatron fluxes contain both primary and 
atmospheric contributions. In order to estimate the primary components 
it is necessary to examine the fluxes as a function of atmospheric depth. 
We have fitted the measured fluxes to functions of the form 

J(T,X) - a(T) p(T,X) + b(T) s(T,X) 
where p and s are the primary and secondary fluxes as a function of 
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atmospheric depth (X) and energy (T) . We have used for s(T,X) the calcu- 
lation by Daniel and Stephens (1974) for zero cutoff rigidity and solar 
minimum conditions. 

The function p(T,X) must be obtained iteratively. A primary 
spectrum P Q (T,0) is assumed (actually, a variety of initial spectra 
were tried) and numerically propagated through the atmosphere using 
bremsutrahlung and ionization energy losses to obtain Pq(T,X). This 
function is then used in a least squares fit to the equation, yielding 
the multiplier a(T) for each ene gy separately. The next estimate 
for the primary spectrum above the atmosphere is thus P^(T,0) • a(T) 
P q (T,X), wnlch can again be propagated through the atmosphere to give 
P^ (T, X) . The P 2 (T,0) obtained from the third least-squares fit to the pre 
ceedlng equation does not differ appreciably from P^(T,0) and is thus the 
best estimate of the primary flux at energy T. Ascent data for several re 
presentative rigidity intervals are shown in Figure 3, along with the 
fitted primary, secondary and total curves. 

In some cases the procedure outlined above gives primary fluxes 
within one standard deviation of zero, some of which are negative (due 
to use of Gaussian statistics). In these cases, the iterative procedure 
has been redone using a Poisson analysis which prohibits negative flux 
values. The results from this calculation are 95% confidence upper 
limits for the desired fluxes. 

IV. RESULTS AND DISCUSSION 

Figure 4 rhows the positron and negatron spectra obtained from 
the analysis described above, and Figure 5 compares our total electron 
spectrum with those from other measurements near solar minimum. 
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The present fluxes fall aubstant lal ly below thoae from Che nearly iden- 
tical experiment carried out In 1972 (Daugherty, 1974; Daugherty et al., 
1975). This might be explained by the lower aolar activity in 1972, as 
indicated by high-latitude neutron monitor levels about 2Z higher than 
in 1974. However, low energy protons observed by the Pioneer spacecraft 
network do not show such an effect. Below about 200 MeV, the present 
psectrum also falls below the 1972 spectrum of Webber et al.(1973), but 
agrees reasonably well with the OGO-5 results of Fulks and Meyer (1973) 
between 50 jnd 200 MeV. 

Fig-ri 6 compares our positron results with those of earlier ex- 
periments. Here again there appears to be some discrepancy with our 
earlier results at about 150 MeV. We note that all of the earlier re- 
sults, including those of Daugherty et al. (1975), used less detailed, 
and therefore presumably less accurate, models for the atmospheric 
electron production and propagation. For this reason, we feel that the 
results of Beuermann et al. (1970), should be considered as upper limits 
to the fluxes between 10 and 50 MeV. In any case, their results are 
consistent with our upper limits, despite having been measured at 
a time approaching the last solar maximum. 

We now consider the implications of these results for the problem 
of solar modulation. The first requirement is a theoretical formulation 
describing the process of solar modulation. The modulation of galactic 
cosmic rays in the interplanetary medium has been discussed in terms of 
a spherically symmetry solar cavity within which the particles undergo 
convection, diffusion, and energy changed as a result of interactions with 
the expanding solar wind and associated magnetic fields. In such a model, 
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Che particle number denelty la determined, in principle, by a Fokker- 
Planck equation where the aolar wind speed, the diffusion coefficient, 
and Interstellar particle spectrum are si'.clfled (e.g., Parker, 196S, 

1966; Cleeson and Axford, 1967, 1968 a, b) . The Integration of the 
Fokker-Planck equation has been achieved for this analysis through a 
numerical technique due to Fisk (1971). A t rlal-and-error method was 
utilized to obtain the modulated spectra for comparison with our mea- 
surements. The solar wind speed was fixed at 400 km/sec and the radial 
distance to the bo idary of the modulating region was taken as 10 AU. 

The radial dependence of the diffusion coefficient, K, was fixed; thus 
the parametric adjustments to find modulated spectra that fit the measure- 
ments involved varying the rigidity dependenc. 1 and multiplicative con- 
stant of the diffusion coefficient. Since the total modulation depends 
principally on the Integral of the diffusion coefficient over the volume 
of the modulating region, the selection of radial dependence and distance 
to the boundary are reasonably arbitrary, and affect only the multi- 
plicative constant. The Interstellar positron and total electron spectra 
are imposed as boundary conditions In order to solve for the modulated 
spectra. 

The Interstellar positron spectrum from Ramaty (1974)was used. 

In the 15 to 800 MeV energy Interval, the positrons result principally 
from collisions of high energy cosmic rays with the Interstellar medium. 

As the high energy cosmic rays are not greatly affected by solar modula- 
tion, their spectra can be accurately decermlned from measurements at 
earth. Thus, the greatest uncertainty in the calculation of the inter- 
stellar positron spectrum, in the energy region of Interest, is the assumed 
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average cosml ray path length, which Ramaty takes as 5 g cm 

The Interstellar electron spectrun le obtained from an analyela 
by Cummings et al. (1973), of measurements of the non-thermal radio 
spectrun. In the direction of the galactic anti-center. The uncertainty 
Is estimated by Cummings et al. , to be a factor of 4 for elect on energies 
above 300 MeV due mainly to uncertainties In the galactic magi itlc field 
strength and the total line of sight emission length. The un ertalntv 
in the electron spectrum Increases toward lower energies, e» eedlng a 
factor of 50 near 100 MeV, primarily due to uncertainties M galactic 
parameters affecting Interstellar radio absorption. 

The positron data are able to most directly addre >s the question of 
solar modulation In this energy region for two reasors: First, the cal- 

culated positron spectrum Is accurately determined as discussed above. 
Second, the positrons In the 20 to 800 MeV energy range measured at 1 AU 
are very probably modulated galactic particles and not of solar or planet- 
ary origin. 

Initial efforts to fit the positron data with a modulated spectrum 
yielded the best fit with a diffusion coefficient proportion to rigidity 
squared. Such a strong rigidity dependence, however. Is In apparent dis- 
agreement with the rigidity dependence as deduced from the locally measured 
magnetic field power spectra (e.g., Joklpll, 1966). 

Further efforts to Improve the fit of the modulated positron 
spectrum to the measurements Included exercising a solar modulation 
model due to Fisk (1975) which utilized off-ecliptlc effects. The modu- 
lated spectrum could then be made to fall more steeply toward lower 
energies while maintaining a diffusion coefficient proportional to 
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rigidity by Invoking a latitude dependent diffusion coefficient and 

solar wind speed. "Reasonable" choices of latitude dependence, however, 

did not result In a dramatic Improvement of the fit of the modul.'ed 

spectrum to the positron data, over the fit obtained with the spherically 

summetrlc model. The spherically symmetric model was then used with the 

constraint that the calculated modulated spectra would provide the best 

mutual fit to both the positron and total electron measurements. The 

results are shown In Figure 7. The interstellar positron spectrum has 

been raised by a factor of 2 which would correspond to an interstellar 

_2 

path length for the cosmic rays of abou*- 10 g cm . With thi'j adjustmen 
to the positrons, the positron and total electron modulated spectra are 
'ti to be reasonably well fit by the measurements above 180 MeV by a 
model using a diffusion coefficient proportional to rigidity, as shown 
by the solid curve in Figure 7. 

At energies below 180 MeV, two approaches are examined. The diffu- 
sion coefficient Is first allowed to retain the same rigidity propor- 
tional form toward lower energy. This yields a modulated positron 
spectrum consistent with the upper limits in the low energy positron 
data; however, the total electron measurements are well above the result- 
ing modulated total electron spectrum. These electrons must then be 
accounted for by a local source. A possible candidate is the Jovian 
magnetosphere, as suggested by Teegarden et al. (1974), to explain 
lower energy quiet time electron increases. The spectral Index of 1.5 
for electrons below 10 MeV, observed by TeeRarden et al., In inter- 
planetary electron Increases, is consistent with the spectral Index 
of the electron excess in our observations. In fact, the two balloon 
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flights which produced the data under discussion in this paper occurred 
just on either side of a large electron increase ns observed by the IMT-7 
low energv electron telescope. 

Historically, ir the absence of a viable local source for nega- 
trons and less accurate galactic positron measurements, the diffusion 
coefficient was postulated to change its rigidity dependence. This 
results in a low energy modulated spectrum as Indicated in Figure 1 
by the dashed line. On the basis of the present data, this possibility 
cannot be excluded. Note that if the transition from a diffusion coef- 
ficient proportional to rigidity, to a constant value, were made in a 
less abrupt manner, the "valley" at about 180 MeV would diminish. This 
would improve the fit to the total electron spectrum but would have 
problems with the 95% confidence upper limit in the positron data at 
150 MeV. 

A more complete understanding of the interplanetary (ar.d galactic) 
posittvu spectra await progress on several fronts. The critical para- 
meters in solar modulation theory have yet to be directly measured out 
of the ecliptic plane, or at radial distances near the bomdary region. 

In addition the plausibility of Jovian electrons with the required 
fluxes above 10 MeV remains unclear. Positron and negatron measurements 
in the 10 MeV to 1 GeV energy interval without the encumbrance of atmos- 
pheric secondaries could eliminate the ambiguity of interpretation still 
possible in this region. 



FIGURE CAPTIONS 


Figure 1 - A schematic view of the Instrument. 

Figure 2 - Total electron fluxes from both flights for final portion 
of ascent and entire float period. 

Figure 3 - Positron and negatron ascent data for representative energy 
Intervals, from both flights combined. 

Figure 4 - PoBltron and negatron spectra obtained from this experiment. 

Figure 5 - Comparison of total electron spectra observed during the 
current solar minimum. 

Figure 6 - Comparison of measured positron srectra betveen 10 and 2000 MeV. 

Figure 7 - Modulation of interstellar spectra of positrons and total 

electrons. Solid line marked "modulated" shows effect of a 
diffusion coefficient proportional to rigidity. Dashed curve 
Indicates the effect of setting the diffusion coefficient 
to a constant value below 180 MeV. Positron and negatron 
data points are from this experiment. 
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